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Abstract 

We propose a continuum model to describe the molecular alignment in thin nematic 
shells. By contrast with previous accounts, the two-dimensional free energy, aimed at 
describing the physics of thin films of nematics deposited on curved substrates, is not 
postulated but it is deduced from the conventional three-dimensional theories of nematic 
liquid crystals. Both the director and the order-tensor theories are taken into account. 
The so-obtained surface energies exhibit extra terms compared to earlier models. These 
terms reflect the coupling of the geometry of the shell with the nematic order parameters. 
As expected, the shape of the shell plays a key role in the equilibrium configurations of 
nematics coating it. 

PACS: 61.30.Dk, 68.35.Md, 61.30.Jf 

keywords: Nematic shells, nematic membranes, two-dimensional nematic order 



1 Introduction 

Nematic liquid crystals are aKg regates of rod-like molecules that tend to align parallel to each 
other along a given direction (|de Gennes and Prosti Il995l ). Due to their easy response to 
externally applied electric, magnetic, optical and surface fields, liquid crystals are of greatest 
potential for scientific and technological applications. Currently, there is an increasing interest 
in soft matter physics on s mall spherical colloidal pa rticles or droplets coated with a thin layer 
of nematic liquid cr ystal (iLopez-Leon et ah . 2011bl lah. The hope is to build mesoatoms with 
controllable valence ( Nelsonl . 2002 ). We refer to these coating layers as nematic shells. 

When nematic liquid crystals are constrained to a curved surface, the geometry induces 
a distortion in the molecular orientation. The possibility to have an in -plane order rather 
than a spatial distribution of the molecules, depends on the sh ell thickness (jVitelli and Nelson , 



20061 : Fernandez- Nieves et al. . 2007 : Lopez-Leon et al. . 2011b ). In ultra-thin shells, the inter- 
action with the colloid surface enforces a sort of degenerate anchoring, i.e. the tendency of the 
molecules to align along any direction parallel to the surface. Thus, unavoidably defects arise 
when nematic order is established on a surface with the topology of the sphere. The number 
of defects is a consequence of the Poincare-Hopf theorem, that states that any configuration 
must have a total topological charge equal to the Euler-Poincare characteristic of the surface. 
For instance on a sphere, whose characteristic is -1-2, we can have two dia metrically opposite 
-1-1 defects or four +1/2 defects located at the vertices of a tetrahedron (|Vitelli and Nelson . 
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20061 ). This tetrahedral defect structure is of gre at interest in material science because defects 
regions can be functionalized to serve as bonds (jNelsonl . 120021) . This could lead to tetravalent 
mesoatoms with sp'^-like directional bonding like carbon. Theoretical studies have emphasized 
the possibility to control t he location of the defects, and hence the valence of mesoatoms, by 
varying the shell geomet r y ( Krali et al.l 120111 ) or by tuning the elastic constants of the nematic 
()Bated . l2008l : IShin et al.l . l2008l) 

Most theoretical studi es on nematic liquid crystal s are framed within the classical director 
theory (see, for instance, (|de Gennes and Prostl Il995t IVirgal Il994l) ). In this setting, the local 
properties of the nematic liquid crystal are described through a unit vector, the director, parallel 
to the local molecular direction. The equilibrium configurations of the nematics minimize the 
Frank's free energy, with respect to all configurations that satisfy the boundary conditions. 
However, the director description of a nematic configuration misses a relevant information at 
the mesoscopic level: the dispersion o f the molecules arouiid the d irector . The order-tensor 
theory, put forward by de Gennes (see de Gennes and Prost 1 1995f ): Virga ([l99i)), overcomes 
this gap by introducing a richer kinematic description. Within this theory the free energy to 
minimize is the Landau-de Gennes free energy. 

Theories for two-dimensional nematic order have been proposed in both director and order- 

tensor schemes (IStralevl. 1971 : Helfrich and Frost . 1988 : Lubenskv and MacKintoshlll993t Selinger et al 



2001 ; Biscari and Terentievl . 2006 ; Krali et al.l . 2011 ) and use free energies derived by symmetry 
arguments or mesoscopic properties. By contrast, our approach derives the surface free energy 
for thin films as limiting cases of the well established three-dimensional theories of nematic liq- 
uid crystals. The main concern is how classical theories (Frank and Landau-de Gennes theories) 
reduce when the nematic molecules, confined within a thin region, align in the direction parallel 
to the underlying surface. A prominent role is played by the ratio between the thickness of the 
shell, denoted by h, and the minimum radius of curvature of the entire shell, denoted by £. In 
fact, the surface versions of Frank and Landau-de Gennes free energies can be deduced from 
the three-dimensional models under the assumption of smallness of the ratio h/£. 

Conversely to existing models, we find that in the two-dimensional directory theory the 
twist term does not vanish. Actually, it expresses the tendency of the molecular axis to 
align with the principal directions. Moreover, our analysis provides a coherent way to ob- 
tain the two-dimensional order-tensor theory. As a result, we retrieve the quadrupolar cou- 
pling between the two-dime nsional order tensor and the curvature tensor already obtained by 
Biscari and Terentiev (l2006l) using mesoscopic arguments. 

The paper is organized as follows. In Section 2, we introduce the mathematical notation 
and terminology. Sections 3 and 4, are devoted to obtain surface free energies from Frank and 
Landau-de Gennes theories, respectively. Mathematical topics employed in these sections and 
some details of the calculations are reported in the Appendixes. Finally, we draw our conclusive 
remarks in Section 5. 



2 Geometrical preliminaries 

We first introduce the terminology and establish some preliminary notations. First of all, three- 
dimensional vectors arc denoted by lower-case boldface letters, whereas second order tensors 
are denoted by upper-case boldface letters. The scalar, vector and tensor products between two 
vectors u and v arc indicated by u • v, u x v and u (g) v, respectively. In cartesian components, 
u-v = UiVi, (ux v)i = eijkUjVk, (u(8)v)y — UiVj, where summation is understood over repeated 
indices, and the third order tensor e = e^k is the Ricci alternator. The composition between 
two second-order tensors A and B is the tensor C = AB with components Cij = AihBhj, 
whereas the composition between a second order tensor A and a vector u gives the vector 
V = Au with components vi = AijUj. Finally, the scalar product between A and B is the 
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Let us assume that the nematic shell occupies a thin region V of thickness h around a regular 
compact surface S. Let Us be the normal unit vector field to S. We parameterize points in the 
bulk through a coordinate set {u, v, ^) such that 

p{u,v,(^)^Ps{u,v)+£,Us{u,v), (1) 

where ps is the normal projection oi p onto S, and |^|, with ^ e is the distance of 

p from the same surface. Such a coordinate set is well defined in a finite neighborhood of S. 
More precisely, we introduce the principal curvatures cis{ps) and C2s{ps) of S at point ps, and 
assume 

/i < min (max{|cis(p5)|, |c2,sb5)|})"^ = ^- (2) 

For every fixed ^ G [— /i/2,/i/2], equation ([!]) defines a parallel surface = {ps + ^i^s{ps) ■ 
Ps & S} at distance |^| from S with the vector field v> : p £ i-^ i^siPs) as unit normal vector 
field. In such a way, the unit vector field is defined on the entire region V. The second-order 
tensor Vz^ is symmetric. Its eigenvectors are f (with a null eigenvalue) and the unit vector 
fields 

ei{p) = ei,{ps) (i = l,2), 

where eis and e2s represent the tangent principal directions fields on S. The spatial gradient 
for each eigenvector is 

Vi^ = — ei(g)ei-- — 62(8)62, (3) 

1 - ^Cis 1 - ^C2s 

Vei = Ki(^)e2 ® ei + ^2(0^2 8) 62 + — -^7 — v ® ei, (4) 

1 - t.Cls 

Ve2 = -Ki(^)ei (8 ei - K2(C)ei (g) 62 + — ^^^r — <^ 62, (5) 
where the functions ki and K2 are given in Appendix A. We refer the reader to the book of 



do Carmd ()1976f ) for a more comprehensive treatise of the geometry of surfaces. 



Let <I> be a smooth field defined on S. Assume $ scalar, vect or or tensor valued. Then the 



surface gradient of $ is defined (see lGurtin and MurdochI (|l975l )) as 



Vs$ = (V$)P, 

where P = I — 14 (8 Mj represents the projection onto the tangent plane of S. The trace of the 
surface gradient of a vector field u defines the surface divergence of u: div^u = trX^u = V,u • P, 
that is a scalar field. On the other hand, the surface curl of u is defined as twice the skew- 
simmetric part of the surface gradient: 

curlgU — — eX^u, 

where e denotes the Ricci alternator. 

Let n denote a unit vector field defined on V such that n(p) = 11(^5) and n- v> = Q at each 
point p. Next, by introducing the conormal unit vector field t = x n, we write the spatial 
gradient of n (see Appendix A for calculation details) : 

Vn = ''^^1 ['^"s ^ ' curls(Lns)] t (g) n + [nt^ — £,Us ■ curls(Lts)]t (g) t 

+ (cn, -C^)i^g'n-r„^i/g)t|, (6) 
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where and tg represent the restrictions of n and t on 5, respectively. The tensor L = —'VjUj 
represents the extrinsic curvature tensor of S. Its trace gives twice the mean curvature H, 
while its determinant gives the Gaussian curvature K. The quantities Cn^ = • Lris, t„^ = 
— • Lng represent the normal curvature and the geodesic torsion of the flux lines of on 5, 
respectively. The latter is zero whenever ng is a principal direction. The quant ities Hn, and 
Kt, denote the geodesic curvature of the flux lines of xis and tg on S, respectively ( Rossol . 2003 : 
Tu and Ou-Yangj . 120041 ) . The geodesic curvature (respectively, Kt^ ) measures the deviance 
of the flux lines of ris (respectively, t^) from following a geodesic on S. Finally, we have set 
i = 1 - 2H^ + Kn^ xi^. 

The divergence and the curl of n are the trace of Vn and the axial-vector associated with 
twice the skew-symmetric part of Vn, respectively. Thus, from © it follows that 



divn 



[nt, - i,i^s ■ curl,,(Lts)] 



(7) 



curln = i ^{ - n - (c„^ - ^K) t + [k„^ - ^z^ • curls(Lns)]i/}. 



(8) 



We observe that the normal curvatures, the geodesic torsion, the geodesic curvatures and the 
surface gradients introduced above are quantities related to the surface S and, therefore, they 
do not depend on ^. Instead, although n has been supposed constant along normal directions 
within the thickness, its spatial gradient depends on ^. 

Finally, since n-^^ = tg ■ (Vns)ns and Kt^ = tg • (V;ns)ts (see Rossol ( 20031 )). the surface 
gradient of is given by 



and consequently 



taVS ® Hs + Kt.tfi 



div,n. 



curLn, 



(9) 



Unlike flat surfaces, the surface curl of ns possesses nonvanishing in- plane components. 



3 Two-dimensional director theory 

The classical elastic continuum theory is based on the pioneering works of Oseen, Zocher and 
Frank p ublished betw een the thirties and the fifties of last century. We refer the reader to the 
book of IVirgal (|1994[ ) for a detailed mathematical treatment. The average alignment of the 



molecules is represented by a unit vector n, called the director, where n is physically equivalent 
to — n. The expression for the elastic energy density (per unit of volume) associated with the 
director distortion consists of four terms 

2wozF = /\i(divn)^ + A'2(n • curln)^ + A'sin x curing 

+ {K2 + is:24)div[(Vn)n - (divn)n] (10) 

where the constants Ki, K2, K3, and K24 are called the splay, twist, bend, and saddle-splay 
moduli, respectively. To ensure a stable undistorted configuration of a nematic liquid crystal 
in the absence of external fields or confinements, the three moduli Ki{i = 1,2,3) must be 
non-negative, whereas the elastic saddle-splay constant must obey Ericksen's inequalities [6]: 

\K24\<K2, K2+K24<2Ki. 

In the absence of external actions, the equilibrium configurations are stationary points of the 
total energy 

W = [ wozF{n,Vn)dV, (11) 
Jv 
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according to the boundary conditions. These may consist in fixing n at the boundary (strong 
boundary conditions) or in allowing n to rotate freely (free boundary conditions) . Intermediate 
situations, known as weak anchoring boundary conditions, can be envisaged by including an 
anchoring energy that penalizes the deviation of the molecules at the boundary from a given 
direction. Furthermore, the free energy density may account for extra terms in order to describe, 
for instance, the interaction of the nematic with external electric or magnetic fields. 

Let us introduce the small parameter s ~ h/£. The smallness of e on the one hand en- 
sures that the parameterization ([T]) is properly defined and on the other hand, with the aid of 
Proposition [2] (see Appendix B), it allows us to approximate the OZF free energy as follows 

Wozf^wSzf^t; [ [fci(div,,n,,)2 + fc2(n,. •curUn,)2 + fc3|n^ xcurUn^PjdA, (12) 

where ki = hKi {i = 1,2,3). Observe that the saddle-splay term has disappeared in this 
approximation since n has assumed to be constant throughout the thickness. In fact, from ([6]) 
it follows 

div[(Vn)n - (divn)n] = tr(Vn)2 - (trVn)^ = 0. (13) 

Comparing equations ([T2|) and ([TO)) we remark that: (i) Wq^p involves a surface integral 
rather than a volume integral, thus we can refer to Wq^p as a surface free energy; (ii) the 
surface elastic constants ki are obtained by multiplying Ki and the thickness h, and, hence, by 
virtue of Ericksen's inequalities, they must be non negative; (iii) the surface free energy involves 
surface differential operators instead of spatial ones. 

It is worth mentioning a peculiarity of curved substrates with respect to planar nematics. 
Unlike the planar case, the twist term cannot be a priori neglected. Indeed, as it has been 
already observed, curls is not orthogonal to n,,. In fact, by using formulae ([9]), equations ([T2| 
reduces to 



WSzF 



which shows that the twist free energy density is proportional to t^__ . The latter vanishes if 
and only if the flux lines of lie along principal dir ections. Thus, th e twist free energy can be 
disregarded whenever spherical shells are concerned ( Shin et al. . 2008|) or whenever the director 



lies along meridians or parallels of an axisymmetric shell ( Chen and Kamien . 20091 ). 

In light of (fH)) . we can give the following intuitive interpretation for the shell- nematic 
interaction. The arrangement of the molecules on a surface is the result of the competition 
between the splay and the bend free energies that try to put the flux lines of rig and tg along 
geodesies of 5, and the twist term that tries to align the flux lines of ns with the curvature 
lines of S. Furthermore, the term proportional to the square of the normal curvature, expresses 
the tendency of the flux lines of to align with the principal direction of minimal curvature. 

From equation (|14p it follows that within the one constant approximation {ki = k2 = /ss = 
k), the surface OZF free energy becomes 

WSzF iV.rispdA. (15) 



A key feature of the free energy (fTSl) is that it differs from the one used in earlier works (jStralev 



19711 : IVitelli and Nelsonl . 120061 : iTu and Seifertl . |2007|) . Indeed, by denoting a the angle between 



the principal direction ei^ and ng, equation (|15p reduces to 

W$zF^\ I {Hc.-u\^ + cl^+T^)AA, 
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where u) represents the spin co nnection field ( Nelson and Pehti , 1987 ; Bowick and Giomi , 2009t ) . 
A glance at equation (21a) of lNelson and Pelitil ( IQStI ) shows that the terms proportional to 
and were neglected. Clearly, this mismatch stems from the fact that free energy den- 
sity in (jl5|) is proportional to the square of the surface gradient of rather than propor- 
tional to the square of co vari ant derivative of ns a s it is customary to assume (see for instance 
(|Nelson and Pehtil . Il987t ) or (|Tu and Seifertl . l2007l )). 



4 Two-dimensional order-tensor theory 

The director theory describes only states with a single optical axis. For closed shells whose 
topology is different from that of a torus, the tangent vector field n exhibits singular points, 
i.e. regions where the local orientational order of the nematic is undefined. As a result, the 
shell often incorporates so-called topological defects. These mathematical singularities can be 
avoided by introducing a tensorial-order parameter, that describes defects as those points in 
which the nematic melts into a liquid phase (isotropic states). Hereinafter we illustrate the 
geometrical meaning of that order parameter. 

We now recall the order-tensor theory for the usual three-dimensional nematics. Let us 
suppose that the orientation of a single molecule is represented by a unit vector m. Microscopic 
disorder is taken into account by introducing a probability measure /p : §2 ^ R+, such that 
/p(m) describes the probability that a molecule placed inp is oriented along m. The orientation 
of the molecular axis is described at each point in space by a point of the unit sphere §^ (or 
by a unit vector). Thus, if H is any subset of S^, the probability of finding in p one molecule 
oriented within Q, is given by 

= / /p(m)da, 

JQ. 

where a denotes the area measure on Nematics posses a molecular mirror symmetry, i.e., the 
head and tail of a molecule can be changed without experiencing any change in the probability 
distribution. Thus, the probability measure is even, /p(m) — /p(— m), and the first moment of 
the distribution fp is zero. 

The second moment of the distribution is the variance tensor M = (m ® m) , where the 
brackets denote averaging with respect to fp. By definition, M is unit trace symmetric and 
semidcfinitc positive. The spectral decomposition theorem ensures that M can be put in the 
diagonal form: 

M = AiGi ® Gi + A2e2 ® e2 + X^e^ ® 63, 

and, since the eigenvalues of M sum up to one, its spectrum is bounded by sp(M) C [0, 1]. 

Nematics may exhibit three different states: isotropic, uniaxial, and biaxial. It is customary 
to define these states by using the order tensor Q = M — il. Thus, we can have: 

{i) the eigenvalues of Q are equal, which implies Qiso = 0; in this case we label the nematic 
as isotropic. 

(ii) At least two eigenvalues are equal, the nematic is called uniaxial. Simple algebraic ma- 
nipulation allows us to write: 

Quru = S (g) U - ^ij . 

The scalar parameter s G [—^,1] is called the degree of orientation, while the unit vector 
u is the optical axis. We retrieve the isotropic phase when s = 0, while the perfect 
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alignment of the molecules is obtained for s = 1. The case s = — i represents flat 
isotropic distributions, in the plane orthogonal to u. 

(Hi) When the eigenvalues of the order tensor are all different, the ncmatic is labeled as biaxial. 
Then we can write the general expression for the order tensor 

Qbia = s ® u - + A (e+ (g) e+ - e_ (X) e_) , 

where A denotes the degree of biaxiality and s e [^5,1] as above. The sign of A is 
unessential, since it only involves an exchange between and e_. The degree of biaxiality 
does always satisfy |A| < i(l — s). Even for biaxial nematic, s = — i represents flat (non 
necessarily isotropic) distributions. 

The fr ee energy comprises two terms: the elastic energy and the Landau-de Gennes potential. 



Following iLonga et al.l (|1987|), the most general quadratic elastic energy can be written as 



Wei(VQ,Q)= / [LiQij^kQij.k + L2Qij.kQik.] + L24{Qij.kQjk - QijQjk.k).i\dV, (16) 
Jv 

where Li, L2 and L24 are constants. Here a comma denotes a partial derivative with respect to 
one of the coordinates. This energy expresses the tendency of the molecules to arrange parallel 
one to each other in a homogeneous state. 

The Landau-de Gennes potential, WlcIG, is a temperature-dependent thermodynamic con- 
tribution that takes into account the material te ndency to spontaneou s ly arr ange in ordered or 
disordered states. Its density is of the form (see de Gennes and Prost ( 1995[) ) 



WLdoiQ) = F{A B, C) + |trQ2 - |trQ3 + ^(trQ^)^ (17) 

where A — Aq{T — Tc)/Tc, Aq is a material-dependent positive constant, T is the absolute 
temperature and Tc is a characteristic temperature; B, C are material-dependent positive 
constants and F{A, B, C) is a positive constant that accounts for the free energy of the isotropic 
phase. We observe that F{A, B, C) plays no role in the minimization of the Landau-de Gennes 
energy density and the stationary points of wlcIG correspond to either isotropic tensors or, 
whenever B^ — 24 AC > 0, uniaxial tensors of the form 

Qcr = s (g) u - il 

with 



^_ B + - 24:AC 
AC ' 

and u G . In addition to the supercooling temperature below which the isotropic state loses 
its stability, there are two other characteristic temperatures for wlcIG'- the nematic-isotropic 

transition temperature 1 H — — ) Tc at which the nematic and the isotropic phase have the 

\ 27 AaC J 

same energy, and the superheating temperature ^1 + 24 A "^"^ above which the isotropic 
phase is the unique stationary point of wr, dG. The resu lting seven characteristic temperature 



regimes for WLdG are discussed in detail by Turzil ( 2007 ) 
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4.1 Degenerate states 

The procedure to derive the two-dimensional free energy for nematic shells is performed in 
two subsequent steps: (a) we have to specialize the free energy to describe degenerate planar 
distributions, where the eigenvector of M with null eigenvalue coincides with i>\ then, (b) as for 
the OZF free energy, we approximate the three-dimensional free energy under the assumption 
of smallness of the parameter e. 

To describe a degenerate anchoring throughout the shell, let us suppose the nematic molecules 
are orthogonal to v and m(p(w, u,^)) = m{ps{u,v)). Since at each point the probability to 
find m in the direction u is zero, it follows that Mi/ = 0. This means that no isotropic spatial 
states are allowed. Let us intro duce n and t the e igenvectors of M orthogonal to u. We write 



the variance tensor in the form (IKrali et al.l . 120111) 



M=i(I-i/(g)i^) + A(n(g)n-t(g)t), 

where A € [— 5, We recognize that two kinds of uniaxial states are allowed: (a) A = then 
i> is the optical axis and the molecules are randomly distributed orthogonally io u] (b) A = ±i 
and the molecules are perfectly ordered along a direction orthogonal to u . The latter case 
coincides with the directory theory analyzed in the previous section. Note that the sign of A is 
inessential since, the order tensors associated with negative values of A and director n coincide 
with the order tensors associated with the positive degree of order —A and director t. 

An alternative and equi valent parameterization of the variance tensor is the following 



(jBiscari and Terentievl . l2006l ) 



M = 9(n®n) + i(l-(7)(I-i/®iv), (18) 

with 9 = 2A e [-1, 1]. 

It is worth noting that this parameterization can be also obtained from the three-dimensional 
order parameter M by imposing s = — i, by taking u along the normal surface and by choosing 
n along one of the two tangential eigenvectors of M. 

4.2 Elastic energy 

Now, let us introduce the traceless tensor Q, associate d with M, in the u sual way: Q = M— il. 



This tensor differs from the one of equation (6) in iKrali et al.l (|2011l ). which is indeed the 
traceless tensor obtained by subtracting from M one-half of the projector on S (which is the 
identity on the tangent plane). 

With the aim of adapting the elastic free energy to the case of degenerate states, we replace 
Q by Q. Since Q and M (as well as Q and M) differ up to a constant, we have VQ = VM; 
thus, in the elastic energy, Q can be replaced by M. 

Furthermore, by using the parameterization (jl8p and with the aid of equations ([3][5l), the 
following identities hold: 

'Mij.k'Mzj,k = 2q^ {(divn)^ + \n x curing + (curln • n)^} + ^iVqp 



2t-^(l - q)iH - K„,0[(1 - (l)H + 2gc„, - (1 + g)K„^ 



xi\ 



^L-^il-q^)K, (19) 
fcM.fej- = (M,j- fcMjfe - MyMjfc,fc),« + M,j^jM,k,k, (20) 
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{Mij^kMjk - M,jMjkM),t = 2qVq ■ [(Vn)n - (divn)n] 

1-q 



div 



2l 



[(1 - q)H + 2qcn, - (1 + q)Kn, xi]u 



(21) 



q^ [(divn)^ + |n x curing] + ^[Vg^ 



qVq ■ [(Vn)n — (divn)n] 



+ r^{l -q){H- K^XW - q)H + 2qc„^ - (1 + g)K„^ xi]. 



(22) 



As for the director theory, in order to obtain the elastic surface free energy, we expand the 
volume free energy as a power series in the small parameter e and consider only the leading 
order term. Thus, by means of Proposition [T] in Appendix B, we obtain 



/ hq'^sq- [nryts ~ Kt^-n.s)AA- / ^3(1 - (7^)K„,dA 
Js Js 



dA 



ih + h- ^h)q^Tl dA, 



(23) 



where /i = h{2Li+L2), h — '^(^2 + 2^24), '3 = /i(2Li+i2+-^24)/2. In the next section we show 
that these elastic constants are subject to restrictions in order to guarantee the positiveness of 
the elastic free energy. 

In order to interpret the contributions of the different terms, we first examine the special 
case where the perfect uniaxial nematic order {q = 1 everywhere) is enforced on the entire shell. 
Equation reduces to 

Wiiq = 1) = / [h « + < + cij - {h +I2- 4/3)t2J dA, 
Js 

that represents a Frank-like surface free energy (to be compared with equation ^^). The ratio 
between the twist and the splay constants can be tuned acting on the constant li. In particular, 
when L2 = 0, then 4Z3 = 2li + I2, and we retrieve the one constant approximation of the Frank's 

energy (HHI). 

By denoting ~ q{ns ® n^) + i(l — q)P the restriction of M to S, the following identity 
holds 



(24) 



the right hand side of this identity is the quadrupolar couphng bet ween the curvature tensor and 
the surface order tensor derived in Biscari and Terentiev ( 2006h employing quasi- microscopic 
arguments. When q is different from zero, this term express the tendency of xig to align along 
one of the two principal directions depending on the sign of the mean curvature. 

The energy term proportional to the square of surface gradient of q clearly expresses the 
tendency of the nematic to arrange in states with constant order parameter. It is worth to note 
that, for topological reasons, states with non zero uniform q are not always allowed. This is the 
case of closed surfaces with the topology of the sphere. 



The term proportional to Gaussian curvature K was already obtained in lKrali et al.l (|2011l ). 
It is a constant term only when q is homogeneous on a fixed surface, by virtue of Gauss-Bonnet 
theorem. 
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Concerning the second integral of the right hand side of , we find the following identity 
(see Appendix C) 



/ g(Vsg) • (K„^t^ - Kt,ns)dA = i / q^iV^a - ■ dl + [ q^KdA, 

Js ^ JdS ^ JS 



(25) 



with a and uj as in previous section. Thus, for closed shells, the density free energy density 
associated with this term is even proportional to the Gaussian curvature. 

4.2.1 Restrictions on the elastic coefficients 

The positiveness of the free energy imposes su i table restrictions to the free energy coefficients. 
Following the approach pursued in lKrali et al.l (|201ll ). let us decompose the surface elastic free 
energy density as follows 

(26) 



■^ii = will + ^i'fi2 + w'fi3> 



with 



eil 



2/ 2 I 2 
9 \l^t, + '^n. 



|V.<z| 



+ hq (^n, tavs - Kt^ns) • VsQ, 



(27) 



f/2 = |[(l + '7) 



2(1 - q^y^ct^ + (1 - qfct \ - h(\ - q')cr,^ct. 



(28) 



w 



el3 



[h - ih +l2- 3h)q^] tI 



(29) 



where the identity = c^^Ct^ — has been used. Then, we recognize that w^^-^ = vi • AiVi 



and w^i2 = '^2 • A2V2, with 



Ai = 



h 


I2/2 


\ 


12/2 


h/4. 



















-I2/2 h/A 1 


= (qKn 




ts,qKt^,Vsq ■ ns) 



Ao = 



1 



h 



h - 2h 



4 V '1 - 2/3 h 

V2 := [(l + g)c„,,(l-g)ctj. 
Hence, it can be easily proved that > if and only if 

/i>0, \l2\<h, 0<h<h, ll+l2<4k, 

or, equivalently, 

Li>0, 2Li + L2>0, IL24I < 2Li +L2, 1^2 + 2L24I < 2ii + ^2. 



(30) 



By assuming Li > and introducing the ratios Ai = L2/L1 and A2 = ^24/^1, the admis- 
sible region in the (Ai, A2)-plane in which the surface elastic energy density (|23p is positive 
semidefinite, is sketched in figure [TJ It is worth noting that the domain in which the clastic 
free energy density (|26p is positive se midefinite c o ntaiii s the domain of nonnegativeness of the 
surface energy density introduced in iKrali et al.l (|201lh . This in turn contains the domain of 
nonnegativeness of the elastic energy density ([T6|. 
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Figure 1: We have set Ai = L2/L1, A 2 = L2AIL1 Si i s the region in which the elastic 
energy density in (|T6|) is nonnegative Cs eelLonga et al. ( 1987 )): Si U 5*2 represent the domain in 
which the surface elastic free energy in Krali et al. (|201l[ )is nonnegative; finally 5i U ^2 U 5*3 
is the region where inequalities the ^ hold. A = (-3/2,1/2), B = (-1,1), C = (1,1), 
(6/5,4/5). 



4.3 Landau-de Gennes potential 

Let us consider the Landau-dc Gennes free energy density (jl7|) . where Q = Q. A straightforward 
calculation gives 



tr(Q 



1 



1 1 2 /TT^N 1 

Following the same arguments given in Appendix B, we readily derive the surface Landau-de 
Gennes free energy 



8- 



q^]dA fore<l. 



where 



F{A,B,C) 



A 
12 



1.^ 
c = n — , 

2 



B 
108 



C 
144 



T~T* 



a = ao- 



r; = 1 + 



B 



C 



ao = hAo, 



We then obtain a Landau-de Gennes potential wit h two constants in which th e cubic term 
vanish es. An analogous expression is proposed in Biscari and Terentjev ( 2006h : iKrali et"aL 
(j201ll ). It is worth noting that, homogenous states with q ^ are allowed only on surfaces 
with zero Euler-Poincare chracteristic. In fact, only in this case it is possible to define a critical 
temperature that generally depends on the shell curvature. 



5 Concluding remarks 

We have deduced the two-dimensional versions of Frank and Landau-de Gennes free energies 
needed to treat the equilibrium of thin nematic films, coating curved surfaces. These models 
have been obtained as limiting cases of the respective three-dimensional models. The formalism 
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proposed applies to rigid shells as well as to flexible surfaces with two-dimensional nematic order. 
Obviously, in the latter case additional energy terms are required to describe the elasticity of 
th e shell. The problem of equ ilibrium can be framed in the general variational scheme proposed 
in iNapoli and Vergori (I2OIOI) . However, the resulting equations for this complex problem are 



strongly non linear and demand a numerical treatment. 

Our rigorous procedure predicts the existence of new terms in the free energy, with respect 
to earlier models. The physical interpretation of these extra terms is widely discussed in Sect. 
3 and Sect. 4. The key results of our analysis are as follows: 

(i) In the context of the director theory for curved nematic thin films, the twist free energy 
does not vanish. This free energy, coupled with the term proportional to , expresses 
the tendency of the molecules to align along the principal direction of the surface with 
minimum curvature. Thus, the extrinsic geometry of the shell influence s the molecu- 
lar alignment in agreement with the results announced in iMbanga et al" ( 2011 ). In a 



forthcoming work, we show how the twist term influences the stability of a nematic on a 
toroidal surface. 

(ii) In the context of Landau-de Genncs theory, we establish a coherent framework to develop 
a two-dimensional order-tensor theory. As a result, we obtain the coupling term (j24[) . 



This term has been already proposed in iBiscari and Terentievi (l2006f ) , but it required an 



additional phenomenological constant in the model. By contrast, since we deduce that the 
coefficient of this energy is the Frank's constant fci, no further phe nomenological cons tants 



should be introduced. We notice that, within the model proposed iKrali et al.l (j2011l ). this 
term docs not appear; this implies the counterintuitive fact that the biaxiality axes can 
be interchanged without affecting the energy. 

Our approach offers the two-fold advantage of being based on well-established theories and, at 
same time, to avoid the proliferation of phenomenological coeflacients in the free energy expres- 
sion. Therefore, our models describe in an economical and exhaustive manner the equilibrium 
configuration of in-plane curved nemati cs. Obviously , our p rocedure can be extended to more 



complicated models as that proposed in iLonga et al.l (|1987[ ). 



We believe that the results outlined in this paper are the basis to study the arrangement 
of two-dimensional curved nematics. We envisage a series of future studies to establish the 
influence of external actions (temperature, electric or magnetic fields), of the shell geometry, 
and of the material coefficients on the nematic shell texture. 
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A Derivation of the spatial gradients of i>', ei, e2 and n 

Let if = ip{u, v) an orthogonal parameterization of S such that 

eis and = e2s- 



Then, for any fixed ^ S [— ft-/2, h/2\, ip^ = (p{u, v) + ^^^(m, v) is an orthogonal parameterization 
of 5^ such that 

Ve,u = - £.cis)v',u and ip^.u = {l-^C2s)'P,v 
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As consequences of assumption ^ , {ei,e2}, with 

eib):=^ = ^-ei,(p5) G V (31) 
\V>i,u\ W.u\ 

and 

e2(p):=^ = ^-e2.(p5) Vp G (32) 
is a local orthonormal basis of the space of tangent vectors X{S^), whereas 

t^yp) = -T = i^siPs) Vp e y (33) 

is the unit normal vector field on S^. We now introduce the following quantities 

/? = • = = -IV^.^ • 

3^ = -t^s,v ■ ^i.v = C2s(l - ^C2s)<P,v ■ ^.v = C2s(l - iC2s)G , 

= • VC/u = (1 - Ccis) V,« • V,u = (1 - Ccis)^^^, 

-P? = VC.u • = (1 - Ccis)(l - Cc2s)V,M • V,!) 0, 

= VC.f • = (1 - Cc2s)^¥',D • V.f = (1 - Cc2s)^G, 

where i? — ip^u • ¥'.« ^-i^^d G = <^.t, • 

We first derive the gradient of v. From ((33|) it follows that 

(Vi^)iv=0, (34) 

and, since it is a unit vector field, 

iv-(Viv)ej = Vi = l,2. 

Moreover, for any fixed ^, — Vi>' restricted to the spa ce of tangent vect ors X(5^) represents the 
extrinsic curvature tensor of S^. Therefore, following Ido Carniol ( 1976 ): 

ei • (V//)ei : 
ei • (Viv)e2 : 
62 • (V//)ei : 

62 • (Vl/)e2 : 









E.G^ 


-F| 1 


- ?Cis 


a^Fi - 

E.G^ 


- F^ 




e^F^ - 
E^G^ 


- F^ ' 




hF^- 




C2s 


E^G^ 


-F| 1 


- Cc2s 



(35) 



by which we deduce that ei and 62 arc the tangent principal directions on S^. Finally, (|34p - p5p 
yield ©. 

Let us now calculate Ve; {i ~ 1,2). From (|3T|) . (|32l) and since (i = 1, 2) are unit vector 
fields, we deduce that 

(Ve,)i/ = = (Ve,)^e, Vi = l,2. (36) 
Next, since {ei, 62, v} is a local orthonormal basis 

(VeO^e, = -(Ve,)^e, V^,J = 1,2, z ^ j, (37) 
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and 



u ■ (Vej)ej = -Gj ■ (Vufei = S.^ 



1 - 



V*,.7 = l,2, 



where 5ij denotes the Kronecker symbol. By means of p7l 



62 • (Vei)ei = -ei • (Ve2)ei = ki(^), 
62 • (Vei)e2 = -ei • (Ve2)e2 = K2(C), 



(38) 



(39) 



where K^ (i) and K2 ( £) are the geodesic curvatures of the hnes of curvature on 5^. Hence, by 
following Ido Carmol ( 1976 ) and since the surface gradient of a scalar- valued function / defined 
in a neighborhood of S may be written as 



VJ 



f,lL 

7^' 



e2,s, 



the geodesic curvatures of the lines of curvature on arc found to be 



E, 



2E^,/G^ 



Kls , ^'^Cis • 628 



(40) 



(41) 



and 



«2(0 



G 



l-eC2s (l-eCls)(l-eC2.)' 



G 



2G5/E^ 2{1 ~icis)GVE (I - (cis){l - ^C2s)Ve 

_ K2s ^V,C2s ■ Bis 

^1~^CU (l-eCl«)(l-eC2.)' 



(42) 



where 



E., 



His — 



and K2s ~ 



G,u 



2EVG 2GVE 
are the geodesic curvatures of the lines of curvature on S. Therefore, ([55 )) " P^ give ^ and 
©• 

We arc now in position to derive the gradient of the director field n. Since n is a unit vector 
field that does not vary with ^ and is pointwisc orthogonal to u, we get 



(Vn)i^ = = (Vn)"^!!. 

Next, we introduce the angle a that n form with ei so that we may write 
n — cos aei + sin ae2 , t = v x n = sin aei + cos ae2 

and 



Vn 



sin aei (X) Va + cos aVei + cos ae2 (8) Va + sin aVe2 . 



(43) 



(44) 



(45) 



Since n and ei are constant throughout the thickness, also the scalar field a satisfies the equality 
a{p) = a{ps) for all p £ V. Therefore, in view of (P0|) the spatial gradient of the scalar field a 
is 



Va = 



=ei 



{l-^Cis)VE {1-^C2s)VG 

Vsa ■ eis VsU ■ e2s 

-61 + 62 



62 



1 - S.C1S 



1 - £,C2s 
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Thus, 



and 



V ■ (Vn)n ~ — n • (V/>')n = 



cis cos^ a + C2s sin a ~ CcisC2s 

(i-ecis)(i-ec2,) 



• (Vn)t = -n • (Vi/)t 
t • (Vn)n 



1 - 2^H + ^'^K' 

{c2s — cis) sin a cos a 



(\^a • eis) cos a + (X^a • e2s) sin a + kis cos a + K2s sin a 



t • (Vn)t 



(i-eci.)(i-ec2.) 

Ci^Kis COS a + C2sK2s sin a + V,C2s ■ ei^ sin a ~ VsCu ■ 62^ cos a 

(l-^ci,)(l-ec2.) 
C2s cosaX^a • ei^ + cis sinaX^a • e2s 

{I - ^Cu){l - ic2s) 

X^a • n^. + Kis cosa + K2s sina — ^divs{c2s sinaei^ — cis cosQ;e2s) 
1 - 2H^ + K„, xi^ 

^n, - £,i^s ■ curls (Ln^) 
1 - 2H£, + xi^ 

— (VsQ! • eis) sina + (X^a ■ e2s) cosa — sina + K2s cosa 

(i-ecis)(i-ec2s) 

— cisKis sina + C2sK2s cosa + VsC2,5 • ei,, cosa + X^ci^ • 625- sina 

(l-Ccis)(l-^C2s) 

C2s sin aVsa • ei^ - cis cosaX7<,a • e2s 

(1-^C1.)(1-CC2«) 

X^a • — Kis sina + K2s cosa — ^divs(c2s cosaeis + ci^ sinae2s) 



1 - 2i7^ + 



- i.i's ■ curls (Lts) 



(46) 
(47) 



(48) 



(49) 



where ris = cos aeis +sin ae2s and tg = i<s x ns arc the restrictions of n and t on S, respectively, 
and L cisCis ® eis + C2se2s ® e2s is the extrinsic curvature tensor on S. The quantities 

Cn, = cis cos^ a + C2s sin^ a and ct, = cis sin'^ a + C2s cos^ a 

are the normal curvatures of the flux lines of ns and ts, respectively, whereas 

Tn, = (cis — C2s) sin a COS a 

is the geodesic torsion of the flux lines of ns . In deriving (|48p and (|49p we have made use of the 
Lioville's formula (see do Carniol ( 1976 ) page 253) for the calculation of the geodesic curvatures 
Kn^ and Ktj, , i. e. 

^ris ~ V,a • ns + Kis cos a + K2s sin a, Kt^, = X^a • tg — kis sin a + K2s cos a, 
and have employed the identity 

divs(i's X u) = —h>s ■ curlsU (50) 
that holds true for any smooth held u defined on S. We may then conclude that yield 
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B Derivation of W^^^ and VFJ^ 

In this section wc shall derive the approximations of the energies pH)) and ([T5| that are valid 
for a homogeneous nematic whenever e ^ 1. 

Proposition 1 Let n and q be smooth fields defined on V . Assume n to be a unit vector field 
such that 

n(p)-jy(p) = ypeV 

and 

n[ps + ^i^siPs)] = n{ps) yps eSy^e [-h/2, h/2l (51) 
and q a scalar-valued field such that 

q[ps + ii^siPs)] = q{ps) e [-1,1] yps e 5,Ve e [-h/2,h/2]. (52) 

Then, denoting by vol{V) the volume of V , 

I 'Mij,kMij^k _ f [(divsiis)^ + |n, x curl,n,p + (cmi^iis ■ n^)^] 



^^^Jv vol(F) Jg area(5) 

2(1 - q)H[{l - q)H + 2gc„J - (1 - q^)K 



2 



-AA, (53) 



5 area(5) 

I M.,-.-M .,..^^ ^ I q- [(div.n.)^ + |n xeurl.n.P]^^ 
e^ojy Yo\{V) Js area(5) ^ ' 



q^sq ■ [(curl^n^ • Us)t - (divsns)!!^] + (1 - q)H[{l - q)H + 2qc„^ 



4 



5 area(5) 

{AUj^kMjk - MijMjk,k),i f 2qV,g • [(curlsn^ • v>s)ts - (divsn^)]!^] 



hm / ^ i — '—^dV = / ^ — -dA 

e^ojy vol(y) ^5 area(5) 

" '^'^-Kn, dA, (55) 



5 2area(5) 



where ^ 

My = griirij + -(1 - q){6ij - ViVj). 



We observe that 

vol(y) = h 
With the aid of equation ([7]) 



area(5) H / dA 

^'^ Js " 



(56) 



,^(div„)^_^^^ / ^^K^f^-ULt.)]-^^ 



VOl(F) y_„/2 ^5, ^'VOI(F) 

de / g^K-e^..-curU(Lt.)]^^^ 

-/V2 ^5 iVol(F) 

- .M^.-^.. -curMLt.)]- 1 

h/2 ^VOl(F) J 
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Since q, Kt^ and Vs ■ curls (Lt^) do not depend on ^, by means of (j56p we deduce that 



''/^ Lh ^area(5) + Y2 / 



area(5) 



uniformly in 5 as £ — > 0. Therefore, recalling ©i, 



lim 



(/^(divn)^ 
e^'o vol(V^) 

We now use equation ([5]) to obtain 
g^(curln • n)^ 



dV 



g^(divsn.,)^ 
5 area(5) 



d^. 



(57) 



YO\{V) 



dV = 



h/2 ^2^2 



iyol{V) 



dA 



and 



g^|n X curlnj' 
vol(F) 



■dV^ = 



/■ f /"''^^ g'^lcn, - + [Kn^ - ■ curls(Lns)]^ \ 



dA. 



Considering that c„^, Tn^, and i/j • curls(Lns) do not depend on ^, by means of (|56p . we 
have 



'^^'^ ih ( area(5) + ^2 ^ '^"^ 



■d^ 



area(5) 



uniformly in 5 as £ — )> 



and 



h/2 



2 i.^ns ■^*) 



"''/^ i/i ^area(5) + J ^^n^dA 

uniformly in 5 as £ — > 0. Thus, from ^ we deduce that 



e//s- curls (Ln,)]2 



area(5) 



and 



lim 

E-i-O 



lim 



q^{n ■ cur In) ^ 
■ vol(F) 

g^|n X curing 
vol(F) 



■dV ^ 



dV = 



g^(ns ■ curls lis )^ 
5 arca(5) 

I lis X curls nsp 
area(5) 



dA 



dA. 



(58) 



(59) 



By following the same arguments which lead to ([57|) - ([5^ and by taking into account that 



Vsg • eis Vsg • e2s 
Vg = — ei + — 62, 

1 - ^Cis 1 - CC2s 



one can easily prove that 



lim 



|V<?| = 



-dV^ 



Hq\- 



lim 



Jv vol(F) area(5) 

{l-q){H-KnM 



dA, 



(60) 



Ojy L^VOl{V) 

'g area(5) 



[(1 - q)H + 2qc„, - (1 + q)K,, 



{l-q)H + 2qCnAdA, 



(61) 
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lim 



[l-q^)K 



dV 



lim 

e-J-O 



gVg • [(Vn)n — (divn)n] 



{l-q^)K 
g arca(5) 



(62) 



dV 



vol(F) 

gVsg ■ [(curlsiis ■ Vs)ts - (divsn^)ns] 
area(5) 



dA. 



(63) 



Now, let us assume now that 5 is a regular surface whose boundary dS is a regular curve, 
and let r be the tangent unit vector field to dS. Then, the normal unit vector field to the 
surface 

Si {pas + ^vsivas) ■ Pas e 95,^ e [-h/2, h/2]} 



IS 



(r - CLt) X u 



Therefore, by means of the divergence theorem we deduce that 
J div|/rl(l - q)[{l - q)H + 2qcn, - (1 + g)K„^ a;i]l/}dF 



l~q 



{l-q)H + 2qcn, -(l + g)K„,- 



(1 - q)H + 2qcn, + (1 + q)Kn, - 



dA 

dA 



/ - q)[{l - q)H + 2qCn^ - (1 + (?)k„^ xi]v ■ NdA 

J Si 

(l-q) (1-(?)F + 2(?c„, -(! + (/)«;„, 
(1 -q) (1 - q)H + 2qc„^ + (1 + q)K, 
By means of (|56p we may conclude that 



dA 

dA : 



h{l - q^)Kn, dA. 



lim / d:w{^:^—$^[{l-q)H + 2qc^,-{l + q)K^^xi]v\dV 



lim 



{l-q^)K 



area(5) + 
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-dA : 



K„ dA 



(1 - q^)K 
s area(5) 



dA. 



We arrive at ([M]) also whenever 5 is a geometrically closed surface, i.e., dS = 0. 
Finally, equations ([SSll-dSn]) immediately follows from and (|57)) - (|M)) . 

From ([55]) it follows that 

lim vol(y) = /iarea(5). 



£^0 



As an immediate consequence of Proposition [T] and (j65p . we have 



M^e/(VQ,Q) = vol(l^)< 



vol(V) 



dV 



vol(F) 



dV^ 



ly,^ for e < 1, 



(64) 



(65) 
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with Wfi as in (j23p . 

On taking </ = 1 in ((57)) - ([59l) we have the fohowing 



Proposition 2 Let n be a smooth unit vector field defined on V such that 



and 
The? 



n{ps + ^i^s) = n{ps) Vp5 e 5, Ve G [-h/2, h/2] 



lim/(^dF=./(^dA (66) 
e^ojy vo\{V) Js area(5) 

^.^ r (n-curln)^ ^ f (n. ■ curl n.)^ 

^^oiv vol(^) 75 area(5) ^ ^ 



|n X curlnp f |n, x curlsiisp 



lim / ^ ^^dy= / ^ TFT^dA. 

=->0Jv vol(V^) Js area(5) 

Therefore, from ([T^ . Proposition [5] and (|55|) . it foUows that 

/• Xi(divn)^ + /V2(n ■ curln)^ + K^ln x curing . 
WozF = vol(y) ^mUy) 

for e < 1, with VK^^^ as in (|T4)) . 



C Geometrical identities 

Let us consider the orthogonal parameterization of S introduced in Appendix A and set 

The metric tensor induced on S by the Euclidean metric tensor, written with respect to the 
system of local coordinates {xi,X2), is 

g = Edx^ (E) dx^ + Gdx^ ® dx^ . 

The Levi-Civita connection associated with the metric g is defined by the Christoffel symbols 



T^2 T^2 AJ? T^l G*,!! G 2 
-•^ 12 - 21 - 777^ - ViiK.2s, 122-^77^- 7^'^2s, i 22 



(69) 

2G 2E /b'""' 2G' 

Then, the (0, 4) curvature tensor of S has components 

^-^ + r,V^.-r^.r^;,j (70) 

= EG(ysK2s ■ eis - VsKls • e2s + njs + '^2s)*^/37<^5p 

= -£'G (i^, • cm\suj)ej3^esp (/3, 7, (5, p,n,X = 1, 2), 
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where e/j^ = 5ipS2-f — (^17(^2^ is the antisymmetric symb ol and a; = — (/tiseis + K 2sS2s) is the 
vector that parameterizes the spin connection Vlp-ys (see iBowick and Giomi that is 



^^,375 = 67 • {Des)ei3 = ujpe^s {13, J, S = 1,2), 



where D ~ is the usual covariant derivative (see iGurtin and Murdoch dQTS )). It is well 



known that the Gaussian curvature of a surface equals the scalar curvature (see do Carniiol 
(|l992h V Therefore 



«„^ = iV^^ = -^fe-curW. (71) 
2^7 

By means of (|7ip we can prove identity (|25p . We first observe that Kn^n^ + Kt,ts = V,a— lj, 
by which Kn^ts — Kt^iis = x (Vja — w), with a as in Section 3. Next, we recall the identity 

IV, .curl,(V3/)=0, (72) 

that is valid for any smooth scalar field / defined on S. Then, applying the surface divergence 
theorem and identities ([SO]) . ([7T|) and ([7^ lead to 

jjli^sq) ■ (Kn,t-Kt,n)dA= i£div,[gV^ x (V,a - a;)]dA 
- i / q^diYslfs X (V^a- u;)]dA 

= ^ / X i^so. - u;)] • kd/ + i / g^w, • curls(V,a - u;)dA 

= i / q^{V,a~u>)-A\+]- [ q^KdA, 
where k is the outward normal to the boundary dS lying on the tangent plane. 
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